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Abstract 

We report the results of a new method of measuring the capillary length for the dendritic crystal growth of non-faceted 
materials. This method uses a nearly spherical crystal held near unstable equilibrium in an oscillating temperature 
field. For the growth of ammonium chloride crystals from aqueous solution, previous published estimates varied by 
over a factor of 20. With this new method, we find a capillary length do of 0.224 + 0.005 nm, similar to that obtained 
for ammonium bromide crystals. 
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1. Introduction 

Dendritic crystal growth is a common form of so- 
lidification observed for the diffusion-limited growth of 
non-faceted materials. It is important technologically 
because many metals and metal alloys solidify in den- 
dritic patterns, but it is also of more general interest as 
one of the standard examples of pattern formation in a 
nonlinear, nonequilibrium system [1]. For reviews, see 
Boettinger et al. [2] and Glicksman and March (3[]. 

In the standard model of diffusion limited growth, 
there are two characteristic length scales [1]. One is the 
diffusion length I, given by I = 2D/v, where D is the rel- 
evant diffusion constant and v is interface velocity. The 
other is the capillary length, do, which is related to the 
energy of the solid-liquid interface. For a pure system, 
do is given by 



do = yT M c p /L 



(1) 



where y is the solid-liquid surface tension, Tm is the 
melting temperature for a flat interface, c p is the specific 
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heat per unit volume, and L is the latent heat per unit 
volume. 

For growth from solution, do is given by 



do = 



7 dfi 
(AC) 2 dC 



(2) 



where AC is the difference in concentration between the 
solid and liquid phases, and pt is the chemical potential. 

For a few simple pure materials, such as helium 
or xenon Q3Q , it is possible to measure the materials con- 
stants in Eq. Q] directly and thus determine do- More 
commonly, however, do must be determined from mea- 
surements of the solid-liquid interface shape under well- 
controlled experimental conditions. 

For pure materials, or for dilute solutions, Schaefer, 
Glicksman, and Ayers J§] showed how to use measure- 
ments of the curvature of the solid-liquid interface near 
a grain boundary to determine do, and applied the tech- 
nique to succinonitrile. Glicksman and Singh f3] ap- 
plied the same technique to pivalic acid. Similar mea- 
surements have been made for xenon fl]] and pyrene 
§]. Luo, Strachan, and Swift [10] developed a method 
to use the maximum supercooling obtainable for a pure 
system to determine the solid-liquid interfacial energy, 
and applied it to the ice-water system. 



For growth from solution, Dougherty and Gollub [ 1 1 1 
modeled the shrinking of a nearly spherical crystal of 
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ammonium bromide. Tanaka and Sano 11211 and Sawada 
et al. [13] applied the same technique to ammonium 
chloride. Gomes, Falcao, and Mesquita [ 14] used the in- 
stability of a cellular interface to measure the capillary 
length for a nematic-isotropic interface in liquid crys- 
tals. 



For ammonium bromide, Dougherty and Gollub 111 111 
estimated do = 0.28nm. For ammonium chloride, 
Tanaka and Sano lfl2ll estimated a much larger value of 
do = 1.59 + 0.06nm, while Sawada et al. [ 13] estimated 
a much smaller value of do - 0.065nm. 



In this paper, we extend the method of Ref. Ill 111 to 
maintain a nearly spherical crystal of ammonium chlo- 
ride close to a state of unstable equilibrium, and extract 
a value for do from the quasistatic growth and shrinking 
of that crystal. Our final result of 0.22nm is intermedi- 
ate between the two previous results, and is also com- 
parable to our previously-reported value for ammonium 
bromide. 



2. Theory 

The theory for diffusion-limited growth of an 
isotropic spherical crystal in an isothermal solution is 
developed in Ref. QJJ] . The central result used here is 
that the radial growth velocity, dR/dt, is given by 



dR _D I 2d 
~dt ~ ~R\ ~ ~R~ 



(3) 



where R is the radius of the crystal, D is the diffusion co- 
efficient for NH4CI in aqueous solution, do is the capil- 
lary length, and A is the dimensionless supersaturation, 
given in Ref. H by 



A = 



ji{T) - ii eq 
AC(dfi/dC) 



In this experiment, the typical temperature variations 
are less than +0.02°C, so we adopt a simple linear 
model for A, namely 



A = c T (T eq - T) 



(4) 



where T eq is the temperature at which a flat interface 
would be in equilibrium, and cj is to be determined em- 
pirically. 

The diffusion coefficient for lower concentrations and 
temperatures was measured by Lutz and Mendenhall 
iflHl and by Hall, Wishaw, and Stokes [ 16]. Extrapo- 
lating those published values to the temperatures and 
concentrations used in this work, we estimate D = 
2500/vm 2 /s, similar to the values of 2600yum 2 /s used by 
Tanaka and Sano 11211 . and 2280jum 2 /s used by Sawada 



et al. 11311 . For this experiment, where the crystal rests 
against the bottom plate of the growth cell, we estimate 
the effective diffusion constant to be 0.71D [ 1 il l 1211 . 



3. Experiments 

The experiments were performed with aqueous solu- 
tions of ammonium chloride with approximately 40% 
NH4CI by weight. The saturation temperature was ap- 
proximately 80°C. The solution was placed in a 40 X 
10x2 mm glass spectrophotometer cell sealed with a 
Teflon stopper held in place by epoxy. A thermistor was 
mounted on the cell for direct readings of the tempera- 
ture near the crystal. The cell was mounted in a mas- 
sive copper block, surrounded by an insulated outer alu- 
minum block, and placed on a microscope enclosed in 
an insulating box. 

The temperature of the outer aluminum block was 
controlled to approximately ±0. 1 °C. The temperature of 
the inner copper block was controlled directly by com- 
puter, allowing complete programmatic control over the 
temperature during the course of a run. The tempera- 
ture of the sample was stable to within approximately 
±5 x 10~ 4 °C. 

Images were obtained from a charged coupled de- 
vice (CCD) camera attached to the microscope and ac- 
quired directly into the computer via a Data Translation 
DT3155 frame grabber with a resolution of 640 x 480 
pixels. The ultimate resolution of the images was 
0.314 ± 0.005/xm/pixel. 

To obtain a single crystal, the solution was heated to 
dissolve all the NH4CI, stirred to eliminate concentra- 
tion gradients, and then cooled to initiate growth. Many 
crystals would nucleate, but an automated process was 
set up to acquire images and slowly adjust the tempera- 
ture until only a single isolated crystal remained. 

The interface position and radius of the nearly spher- 
ical crystals were measured as in Ref. 01711 . The temper- 
ature was then continually adjusted until this crystal was 
the desired size. Although such a spherical crystal is in 
a state of unstable equilibrium, we found it possible to 
maintain it indefinitely, provided we continually moni- 
tored the size and adjusted the temperature accordingly. 
This isolated, nearly spherical, crystal was allowed to 
stabilize for several days, and is shown in Fig.Q] 

Once a single, stable crystal was obtained, the tem- 
perature of the sample was slowly oscillated with an am- 
plitude of 0.0125°C and a period of 1 hour. This caused 
the crystal to alternatively grow and dissolve slowly, 
such that there was plenty of time for the system to at- 
tain thermal equilibrium, but the size and growth rate 
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Figure 2: Measured (dots) and fit (line) values for the radius as a 
function of time of an oscillating crystal. The driving temperature 
variations AT are shown with a dashed line (red online). 



Figure 1: Crystal of ammonium chloride held in unstable equilibrium. 
The image is 100 (im across. 



changed enough that both terms on the right of Eq. [3] 
were significant. 

We then used the data for T(f) to numerically inte- 
grate Eq. [3] and compare to the measured R(t). We de- 
termined the constants cj, T eq , and do in Eqs. [3] and [4] 
by minimizing the squared difference between the 
integrated prediction and the measured values for R(t). 
Since both cj and do are ultimately multiplied by D in 
Eq. [3] these measurements only constrain the products 
Dcj and Ddo- For simplicity, and for consistency with 

Refs. dm m m. we have assumed an extrapolated 

value for D, but this experiment provides no indepen- 
dent confirmation of that value. 



3.1. Results 

The results are shown in Fig. [2] The best fit param- 
eters are d Q = (2.24 + 0.05) x 10~Vm, c T = (4.23 + 
0.01) x 1(T 3 / C, and T eq = 80.70 + 0.01°C. A second 
run on a different crystal gave consistent results. 

This value compares favorably with the value of 2.8 x 
10 _ V m previously reported for NF^Br and is in 
between the values for NH4CI reported in Refs. lfl2"ll and 
10. 

With this value of do, we can also make an updated 
estimate for the stability constant <x*, defined by 

. 2d D 

°" = — r 



Based on the previously reported [18] value of vp 2 — 
12 +2^m 3 /s, we estimate <x* = 0.093 +0.016 for growth 



at small supersaturation. 

There were two main factors we found to be impor- 
tant in getting consistent results from this experiment. 
First, control and measurement of the temperature at the 
sample was critical. In preliminary trials with poorer in- 
sulation, or where we measured the temperature of the 
surrounding copper block instead of the sample cell, we 
were unable to get satisfactory fits to Eq. [3] over long 
time scales as in Fig. [2] 

Second, the oscillating protocol proved a robust 
probe of Eq. [3] In previous work on slowly dissolv- 
ing crystals, as in Ref. OlUl . we typically found that the 
fitted value for do depended very sensitively on the last 
few data points, where the crystal was shrinking rapidly 
just before completely dissolving. Small drifts in tem- 
perature on the order of +0.00 1°C could make signifi- 
cant differences in the quality of the fit. In Ref. 11211 . the 
fit does not include the last two data points, though it is 
not clear why. By contrast, Fig. [2] includes both grow- 
ing and dissolving segments so that both terms in Eq. [3] 
were significant, and the fit was well-constrained. 

One remaining uncertainty is the possible role of con- 
vection in this experiment. The density differences 
caused by the growth and dissolution of the crystal 
could drive small convection currents. However, the ex- 
cellent fit to the data suggests that Eq.[3]still provides an 
effective model for the growth process. 
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